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We explore the analogy between quark confinement and the Meissner effect in 

superconductors. We measure the response of color-magnetic "supercurrents" from 
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0^ . Dirac magnetic monopoles to the presence of a static quark-antiquark pair in four- 
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dimensional U(l) lattice gauge theory. Our results indicate that in the confined phase 
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these currents screen the color-electric flux due to the quarks in an electric analogy 
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of the Meissner effect. We show that U(l) lattice gauge theory obeys both a dual 
London equation and an electric fluxoid quantization condition. 
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Since free quarks have never been isolated, there must be a mechanism for permanently 
confining them within hadrons. This mechanism is expected to arise in a simple way from the 
theory of strong interactions. It was suggested many years ago [[I]] that confinement would 
occur if the vacuum reacted to a color electric field in a manner similar to the response of a 
superconductor to an external magnetic field. Because of the Meissner effect, if two magnetic 
monopoles of opposite magnetic charge are introduced into the interior of a superconductor, 
the Cooper pairs give rise to persistent currents to generate a counter magnetic field to 
expel the magnetic flux. As a result, the magnetic flux lines from one monopole to the other 
are confined to a narrow Abrikosov flux tube which is surrounded by persistent circulating 
currents. The energy of such a configuration is proportional to the separation of the two 
monopoles, thus permanently confining them. By analogy, therefore, if the vacuum naturally 
expels a color electric flux, the field lines from a static quark-antiquark (qq) pair would not 
spread out in a dipole field pattern but would instead form a narrow flux tube, leading to a 
quark potential proportional to their separation and confinement. 

For this mechanism to work the vacuum must contain objects that react to a color 
electric field in a fashion similar to the reaction of the Cooper pairs in a superconductor 
to an ordinary magnetic field. One possibility is to mimic [TJ the Ginzburg-Landau theory 
of superconductivity by adding to the gauge theory an elementary Higgs field to act as the 
superconducting order parameter, but there is no experimental evidence for any elementary 
scalar particles in particle physics. The dual superconductor mechanism |2| is an alternative 
that does not require the ad hoc introduction of a Higgs field but instead uses dynamically 
generated topological excitations to provide the screening super currents. For example, U(l) 
lattice gauge theory contains Dirac magnetic monopoles in addition to photons 0,|]]. The 
dual superconductor hypothesis postulates that these monopoles provide the circulating 
color magnetic currents that constrain the color electric flux lines into narrow flux tubes, 't 
Hooft has shown |J that objects similar to the Dirac monopoles in U(l) gauge theory can 
also be found in non-Abelian SU(N) models. 

U(l) lattice gauge theory in 4 dimensions has both a confined phase at large charge 



and a weak coupling deconfined phase corresponding to continuum electrodynamics with a 
Coulomb interaction between static quarks. Therefore confinement or its absence can be 
studied using U(l) lattice gauge theory as a prototype, before tackling the more complicated 
non-Abelian theories that actually describe quarks. Much evidence for the dual superconduc- 
tor hypothesis has accumulated from studies 0||J6|-0 of U(l) lattice gauge theory. Polyakov 
|| and Banks, Myerson and Kogut Q showed that U(l) lattice gauge theory in the presence 
of a quark- ant iquark pair could be approximately transformed into a model describing mag- 
netic current loops (the monopoles) interacting with the electric current generated by the qq 
pair. DeGrand and Toussaint || demonstrated via a numerical simulation that the vacuum 
of U(l) lattice gauge theory was populated by monopole currents, copious in the confined 
phase and rare in the deconfined phase. This behavior has also been seen in non-Abelian 
models after gauge fixing [HJ . Many studies of non-Abelian models using Dirac monopoles 



10|— L2| or other topological excitations [13| support the dual superconductor mechanism, 



although other studies |b|] dissent. 

So far, studies of confinement have examined "bulk" properties such as the monopole 
density J6| |iTi| , |To| , the monopole susceptibility ]T|,|TT| and the behavior of the static quark 



potential P,|TT||. However, the case for the dual superconductor hypothesis is incomplete 
without an explicit demonstration that a static qq pair actually induces the appropriate 
persistent current distribution. In this Letter we present the first direct evidence for this 
behavior. We further show that there are exact U(l) lattice gauge theory analogues of two 



key relations [IB] that lead to the Meissner effect in a superconductor; the London equation 
and the fluxoid quantization condition. 

Our simulations are done on a Euclidean spacetime grid of volume L 3 x L t , where L is 
the spatial size and L t the temporal size of the lattice in units of the lattice spacing a. The 
U(l) gauge degrees of freedom are complex numbers of unit magnitude residing on the links 
of the lattice and are written U^if) = exp[z# M (r)], where r denotes a point on the lattice 
and \i the direction of the link from that point. The links form a directed lattice so that 



U-^if+n) = UUr). We use a standard Wilson action S supplemented with a Wilson loop 
W to represent a static qq pair with charges ±1. 

s = S - iW 
= f3 £ [l-cosMOME'WW (!) 

f,fl>U ffl 

Here /3 = frc/e 2 is a dimensionless measure of the strength of the charge and exp[i6 , Atjy (r)] = 
U f i(r)U u (f+n)UJ l (f+u)Ul(r) is an oriented product of gauge variables around an elementary 
plaquette of the lattice. The current J^if) is ±1 along the world line of the qq pair and 
otherwise. In the naive continuum limit (a — > 0) S reduces to the action for a pure photon 
field in the presence of a current loop. Physical observables are given by expectation values 

(A) = -Tr (e~ s ° +iW A) (2) 

where Tr denotes an integral over all angles O^r) and Z = Tr exp(— S +iW) is the partition 
function. 

The two observables we study here are the electric flux through a plaquette, which in 
lattice variables is E^if) = Imexp[i^ 4 (r)], and the curl of the monopole current density. 
The latter is found by a prescription devised by DeGrand and Toussaint || , which employs 
a lattice version of Gauss' Law to locate the Dirac string attached to the monopole. The net 
flux into each plaquette is given by {O^if) mod 2tt). If the sum of the fluxes into the faces 
of a 3-volume at fixed time is nonzero, a monopole is located in the box. The net flux into 
the box at fixed time thus yields the monopole "charge" density, or the time component of 
the monopole 4-current Jm, measured in units of the monopole charge gu = 2nhc/e. The 
spatial components are found in a similar manner. The monopole currents form closed loops 
due to the conservation of magnetic charge. 

Our simulations are performed on a 9 3 x 10 lattice using skew-periodic boundary condi- 
tions. Less extensive work on a 7 3 x 8 lattice yields similar results except for the expected 
increase in statistical fluctuations arising from the smaller lattice size. We used a 3 x 3 
Wilson loop oriented in the zt plane and measured the electric flux and the monopole cur- 



rent in the transverse (xy) plane midway along the axis connecting the qq pair. A standard 
Metropolis algorithm |I7J alternated with overrelaxation |18j is used to generate configu- 
rations distributed according to exp(— S ). In the confined phase, we thermalize for 10000 
sweeps and sample the data every 10 sweeps for a total of 7000 measurements, which are 
then binned in groups of 5. In the deconfined phase only half as many measurements are 
taken since the fluctuations are much smaller. Because of the geometrical symmetry of the 
measurements only the ^-components of (S) and (V x Jm) are nonzero. If the Wilson loop 
is removed, even the ^-components average to zero, so the response is clearly induced by the 
presence of the qq pair. 

Figure [3](a) shows the electric flux distribution for (3 = 1.1 where the vacuum is in the 
deconfined phase. The broad flux distribution seen is identical to the dipole field produced 
by placing two classical charges at the quark positions, except that the classical value of 
the flux on the qq axis is a factor of two smaller. We measure the total electric flux from 
one quark to the other, including not only the flux through the plane between the charges 
(0.8504 ±0.0045) but also the flux (0.0951 ±0.0028) that flows through the lattice boundary 
because of the periodic boundary conditions. This yields a total flux of (0.9453 ± 0.0053), 
close to the theoretical value $ e = e/VKc = l/v7? = 0.9534. 

Figure |l|(b) shows the electric flux in the confined phase {(3 = 0.95). In this case the 
flux is confined almost entirely within one lattice spacing of the axis and essentially no flux 
passes the long way around through the lattice boundary. The net flux is again equal to 
l/V/3 within statistical error. This behavior is exactly what one would expect from the 
superconducting analogy, where the flux has been "squeezed" into a narrow tube. The data 
in Fig. [L] are consistent with flux profiles found in other U(l) studies |i| and studies of 
SU(2) Eg and SU(3) 0. 



The electric field E produced by the monopole currents arises from the dual version of 
Ampere's Law 

-cVxE = J M (3) 



By analogy to the superconductor, in order to see a dual Meissner effect we should see a 
dual London relation [|11J between the field and the current of the form 

A 2 

E = — V x J M (4) 

c 

where A is the "London penetration depth" for the electric field. These two equations result 
in electric fields being confined to a region of size A. 

We show in Fig. ||(a) (£) and in Fig. ||(b) — (V x Jm) in the confined phase as a function 
of the distance from the qq axis. The data show that the spatial variation of the flux and the 
curl of the current are very similar, except for the point on the axis which will be discussed 
below. Figure ||(c) shows the best fit found for (£) — (A 2 /c)(V x Jm)- We find a value of 
A/a = 0.482 ± 0.008, which is consistent with the range of penetration of the electric flux in 
Fig. |2](a) and the thickness of the current sheet in Fig. 0(b). The dashed curve in Fig. |2| is the 
result of using the continuum Eqs. (|) and (f|), and it agrees well with the flux distribution 
from the lattice simulations. We also expect that, as in a superconductor, the transition to 
the deconfined phase will be signalled by a divergence of the London penetration depth. We 
have therefore measured A further from the deconfinement transition at f3 = 0.90, and find 
a smaller penetration depth of A/a = 0.32 ±0.02. In the deconfined phase we find an almost 
insignificant value of (V x Jm) and fitted values of A were larger than our lattice size. 

The anomalous behavior of the point on the qq axis can be understood by recalling that 
a superconductor penetrated by an Abrikosov flux tube becomes multiply connected and 
the London relation is replaced by the more general quantization relation for the fluxoid. 
Since our U(l) vacuum is pierced by an electric flux tube, we expect a dual version of the 
fluxoid quantization relation to hold 

JE-dS £j M -d£ = n® e , (5) 

where n is an integer and $ e = 1/V3 * s the quantum of electric flux. In fact, the data in 
Fig. 0(c) represent a lattice version of a delta-function whose strength (1.016 ±0.014) is very 
close to $ e = l/\//3 = 1.026. Thus, if the surface integral in Eq. @ includes the axis of the 



qq pair, we obtain n = 1, while if the axis is excluded from the integral we obtain n = and 
Eq. (^) holds. We examined the response to doubly-charged quarks to find n = 2 electric 
flux quantization but did not get data of sufficient quality to draw any conclusions. 

Equations (^) and (|^) show that, except for the interchange of electric and magnetic 
quantities under duality, the confined phase of U(l) lattice gauge theory behaves exactly 
like a superconductor in an external magnetic field. It is perhaps surprising that a nonlinear, 
strongly interacting, model such as U(l) lattice gauge theory could be described by such 
a simple model as the linear London equations, but our results indicate that the operators 
(£) and (V x Jm) are clear indicators of the presence of confinement of electric flux by 
monopole currents as the dual superconductor mechanism supposes. Although the Meissner 
effect itself requires only that Eq. (|J) hold, our data also support the more restrictive fluxoid 
quantization relation (||). Thus the analogy between Cooper pairs in a superconductor and 
magnetic monopoles in gauge theory is very strong. Analogous studies of SU(2) lattice gauge 
theory along these lines are currently in progress. Because the monopoles appear pointlike 
in our simulations, lattice gauge theory looks like an extreme type-II superconductor and it 
is tempting to argue that the phase transition in lattice gauge theory is a Bose condensation 
of magnetically charged particles, similar to the Bose condensation of charged local pairs 
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FIGURES 
FIG. 1. Surface plot of the electric flux through the xy plane midway between the qq pair when 

the system is in (a) the deconfined phase {/3 = 1.1) and (b) the confined phase (/3 = 0.95). The 

line joining the pair is located at (0,0). 

FIG. 2. Behavior of (a) the electric flux, (b) the curl of the monopole current and (c) the 
fluxoid in the confined phase {(3 = 0.95) as a function of the perpendicular distance R from the qq 
axis. The dashed line in (a) shows the flux expected using the continuum relations (^) and (||). 
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